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Abstract. This paper is concerned with the study of solutions to discrete 
parabolic equations in divergence form with random coefficients, and their 
convergence to solutions of a homogenized equation. In rate of conver- 
gence results in homogcnization and estimates on the difference between the 
averaged Green's function and the homogenized Green's function for random 
environments which satisfy a Poincare inequality were obtained. Here these 
results are extended to certain environments in which correlations can have ar- 
bitrarily small power law decay. Similar results for discrete elliptic equations 
were obtained in 1121 . 

1. Introduction. 

In this paper we continue the study of solutions to divergence form parabolic 
equations with random coefficients begun in [11J . In |1 1 j we were concerned with 
solutions u(x, t, u>) to the equation 

(1.1) M^t,uj) = _ v * a ^ Tx itW ) Vu ( IitiW ) f x e z d , t>o, uj en, 

with initial data 

(1.2) u(x,0,u) = h(x), xeZ d 1 Luen. 

Here Z d is the d dimensional integer lattice and (fi,.F, P) is a probability space 
equipped with measure preserving translation operators T x ,t '■ fi — > 0, x G Z d ,t G 
R. In (jl.ljl we take V to be the discrete gradient operator defined by 

(1.3) V<t>(x) = (Vi#c),... V^(aO), V i <t>(x)=<t>(x + e i )-<f>(x), 

where the vector G Z d has 1 as the ith coordinate and for the other coordinates, 
1 < i < d. Then V is a d dimensional column operator, with adjoint V* which is a 
d dimensional row operator. The function a : — > R d ( d + 1 )/ 2 from f2 to the space 
of symmetric d x d matrices satisfies the quadratic form inequality 

(1.4) \I d < a(w) < AI d , w G n, 

where Id is the identity matrix in d dimensions and A, A are positive constants. 

One expects that if the translation operators t x j are ergodic on O then solutions 
to the random equation ()l.ip converge under diffusive scaling to solutions of a 
constant coefficient homogenized equation. Thus suppose / : R d — > R is a C°° 
function with compact support and for e satisfying < e < 1 set h(x) = f(ex), x G 
Z d , in (|1.2p . and let u e (x,t,uj) denote the corresponding solution to (|1.1[) with this 
initial data. It has been shown in [27j , just assuming ergodicity of the translation 
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operators, that u £ (x/e,t/e 2 ,uj) converges in probability as e — > to a function 
Uh om (x, t), x G R d , t > 0, which is the solution to a constant coefficient parabolic 
PDE 

(L5) du h Ux,t) = _ v * ahomVuhom(Xit) ) j G R d t > 0j 

with initial condition 

(1.6) UhomfoO) = /(a;), x £ R d . 

The d x d symmetric matrix ah om in (|l-5[) satisfies the quadratic form inequality 
(|1.4p . Similar results under various ergodic type assumptions on fl can be found 
in [H [8] HH [35] . In time- independent environments the corresponding results for 
elliptic equations in divergence form have been proven much earlier -see [251 HS1 133 

HE]. 

In [llj we were concerned with obtaining a rate of convergence for the homog- 
enized limit, lim £ _j.o u e (x/e, t/e 2 , uj) = Uhom^j t). The corresponding problem for 
elliptic equations has been extensively studied, beginning with the seminal work of 
Yurinskii [37] ■ Recent papers on the subject have addressed the issue of obtaining 
optimal rates of convergence [HI HOI HE] , and include results for fully non-linear el- 
liptic equations [5] . Optimal estimates on variances of solutions have been obtained, 
but precise results on fluctuations analogous to the central limit theorem have been 
proven only in the case of one dimension [2]. In all these papers one must make 
a quantitative strong mixing assumption on the environment (fi, J 7 , P) in order to 
obtain a rate of convergence in homogenization. For the parabolic problem we were 
unable to find in the literature any results on rate of convergence in homogeniza- 
tion, except for the recent preprint |30j (see also |21| ) in which the environment 
is fixed in time. In [TT] as in [30] our results are restricted to obtaining a rate of 
convergence for the mean (u E (x/s, t/e 2 , ■)) of the solution of (|1.1[) to Uh om (x, t). We 
were able to show that, for certain environments (f2, J 7 , P) satisfying a quantitative 
strong mixing condition, there exists a > depending only on d, A/A such that 

(1.7) sup |( u e (x/e,t/e 2 , ■) ) - u hom (x,t)\ < Ce a for < e < 1. 

xGeZ d J>0 

In [TT] we followed the approach of Naddaf and Spencer [32] to the problem of 
obtaining rates of convergence in homogenization by formulating the quantitative 
strong mixing assumption on the environment as a Poincare inequality. Specifically, 
consider a measure space (f2, J-) of time dependent vector fields u : Z d x R ->• R fc 
with the property that the functions t —± G)(x,t), f 6 R, are continuous for all 
x G 7i d and each oj(x,t) : O — > R fc is Borel measurable with respect to the a- 
algebra T . For a function G : (l — > R the gradient of G is defined in a weak 
sense. Thus if h : Z d x R — > R fc is continuous with compact support the directional 
derivative dhG(Q) of G(Cj) in the direction h is defined as the limit 

(1.8) d h G(u) = lim[G(w + 6h) - G(u)]/S . 

5- fO 

The function da,G(oj) : Z d x R — > R fc is then the gradient of G at w if it is Borel 
measurable and 

/•oo 

(1.9) d h G(u) = dtd Q G(x,t;u>)-h(x,t) = [d Q G(ui),h] 
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for all continuous h : Z d x R — > R fc of compact support. In (|1.9j) we have denoted 
by [•, •] the Euclidean inner product on L 2 (Z d x R, R fc ). Letting || • H2 denote the 
corresponding Euclidean norm, a probability measure P on (O, IF) is said to satisfy 
a Poincare inequality if there is a constant Kp > such that 

(1.10) Var[G(-)] < Kp{ \\dzG(-,Q)\\l ) for all C 1 functions G : Q, -> C. 

If the translation invariant probability measure P is Gaussian, then the measure 
is determined by the 2-point correlation function T : Z d x R — > R fe g) R fe defined 
by T(x,t) = ( Co(x,t)Co (0,0)* ), x G Z d ,t G R, where £(■,■) G R fe is assumed to 
be a column vector and the superscript * denotes adjoint. Defining the Fourier 
transform of a function h : Z d x R — > C by 

/oo 
/l(x, t)e lX<+lte , C 6 [-7T, 7T] rf , G R , 
xe^ -°° 

one can easily see that the Poincare inequality (|1.10[) holds if and only if T G 
L°°([-7T, 7r] d x R). Hence if T is intcgrablc on Z d x R then (fl~T0[) holds, but it is 
unlikely to hold if T is not integrable. 

In the present paper we shall prove rate of convergence results in homogenization 
of the parabolic PDE (jl.l| for certain environments that include some Gaussian 
environments in which T is not integrable. To do this we extend the method 
introduced in [12] for elliptic PDE in divergence form to the parabolic case. The idea 
is to consider environments defined by a(w) = a(cj(0, 0)) where ui : Z d xR4 R™ is a 
translation invariant function of Q : Z d x R — > R k . The gradient of uj with respect 
to Co is assumed to satisfy a uniform integrability condition, and the probability 
space (ft, J 7 , P) for Co to satisfy the Poincare inequality (jl.lOp . 

We define what we mean by the terms used in the previous paragraph. Let 
(il, J 7 , P) be the probability space for uo induced by the probability space J 7 , P) 
for Co and the functional dependence Co —> uo. Translation operators r X:tl x G Z d , t G 
R, on Vt are defined by T x .tto(z, s) = lo(x + z,t + s), z G Z d , s G R, with a similar 
definition of translation on Q. The function Co — > to, which we denote by w(-, - ,Co) 
is translation invariant if 

(1.12) t x jUj(-,-,Co) = lo(-,-,t x jCo) for x G Z d ,t G R, Co G 

Note that if co is a linear translation invariant function of Co then uo is the convolution 
of some function h : Z d x R — >• R" ® R fe from Z d x R to n x k matrices with Co, 
(1.13) 

/oo 
/i(x —y,t— s)Co(y, s), x G Z d , t £ R. 

For given x G Z d ,t G R we use the notation of ([1.9)1 to write the gradient of the 
function oo(x, t, •) : f2 — )• R™ as d^uo(z, s; x, t, Co), z G Z d , s G R, w G f2. The uniform 
integrability condition is then that 

1 

< (i^ 9 ) 9 <co 

for some q with 1 < q < 2. It follows from ([1.13)1 that when w is a linear function 
of Co the condition (|1.14|) is equivalent to the condition that the function h in ([1.13)1 
is q integrable. 



(1.14) 



E 

Tf=Z d 



dt 



sup \dcjLo(0,0;x,t,Co)\ 
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In [IT] we proved that ()1.7|) holds in the case where (fi, J 7 , P) is the stationary 
process associated with a massive Euclidean field theory. This Euclidean field 
theory is determined by a potential V : R d — > R which is a C 2 uniformly convex 
function, and a mass m > 0. Thus the second derivative a(-) = V"(-) of V(-) is 
assumed to satisfy the inequality (jl.4l) . Consider functions cf> : Z d x R — )• R which 
we denote as </>(x, t) where a; lies on the integer lattice Z d and i on the real line R. 
Let f2 be the space of all such functions which have the property that for each x £ Z d 
the function t — > <f>(x,t) on R is continuous, and J- be the Borel algebra generated 
by finite dimensional rectangles {</>(•,■) 6 ^ : \<fi(xi,ti) — Oj| < i = 1,...,JV}, 
where (x,,^) e Z rf x R, a, e R, n > 0, i = 1,...,N, N > 1. For any d > 1 
and m > one can define [H I17j a unique ergodic translation invariant probability 
measure P m on (fi, J 7 ) which depends on the function V and m. In this measure the 
variables <fr(x,t), x £ Z d ,t > 0, conditioned on the variables c/>(x,0), x £ Z d , are 
determined as solutions of the infinite dimensional stochastic differential equation 
(1.15) 

d(j)(x,t) = - . V !-{V(V(t)(x',t))+m 2 <f)(x',t) 2 /2}dt+dB(x,t) , leZ^OO, 
ocbtx.t) £ — ' 2 
v ' x'ez d 

where B(x,-), x £ Z d , are independent copies of Brownian motion. Formally 
the invariant measure for the Markov process (|1.15l) is the Euclidean field theory 
measure 



(1.16) exp 



- ^ V{V4>{x)) +m 2 4>(x) 2 /2 

xEZ d 



[ d<fi(x) /normalization. 

x£Z d 



Hence if the variables <j>(x, 0), x £ Z d , have distribution determined by (|1.16p . then 
</>(-, t), t > 0, is a stationary process and so can be extended to all t £ R to yield 
a measure P m on (57, J 7 ). The probability space (J7,J-", P m ) satisfies the Poincare 
inequality (jl.lOj) with constant Kp m = 4/m 4 . In [TT| we conclude from this that 
the inequality (fi~7|) holds provided a(w) = a(0(0, 0)), where a : R R d ( d + 1 )/ 2 is 
assumed to be a C 1 function satisfying (|1.4[) and ||£)a(-)|| O0 < oo. 

Let (f2, T, P) be a probability space for which the Poincare inequality (|1,10[) 
holds, and uj — > u a function which satisfies the translation invariant condition 
(|1.12[) and the uniform integrability condition (|1.14j) for some q with 1 < q < 2. 
Our goal in the current paper is to show that the inequality (|1.7[) holds for the 
environment (f2, J", P) of w £ fl where a(w) = a(w(0, 0)) and a : R — > R, d ( d +!)/ 2 i s a 
C 1 function satisfying (|1.4j) and HD^^Hoo < oo. Rather than attempt to formulate 
a general theorem for such environments, we shall only rigorously prove that (|1.7[) 
holds for certain limits of the probability spaces (fi, J", P TO ) defined by (|1.15|) . (|1 . 16[) 
as m — 0. In §2 we indicate the generality of our argument by showing that the 
proof of Proposition 6.3 of [11] formally extends to the environment (fi, J 7 , P). 

From (|1.15|) we see that the stationary process u>(-, •) = <p(-, •) is a translation 
invariant function of the white noise stationary process •) = dB '(■,■). It is 
well known that the white noise process satisfies a Poincare inequality (|1.10[) with 
Kp = 1. Consider now the terminal value problem for the backwards in time 
parabolic PDE 

(1.17) = lv*V"(V<f>(z,s))Vu(z >S ), s<t,z£Z d , 
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u(z,t) = uo{z), z G Z d , 

with solution 

(1.18) u(z,s) = G(z, s;x,t, 4>)ug(x) , s <t,z &Z d . 

xGZ d 

We see from pj] that the gradient of uj(x,t) = <j>(x,t), x G Z d ,t G R, with respect 
to Gj should be given by the formula 

(1.19) dzuj(z,s;x,t,u>) = e" m2(t ~ s)/2 G(z, s; x, t, 4>) for s < t, 
cIquj(z, s; x, t, £>) = for s > t, z G Z d . 

In |18j a discrete version of the Aronson inequality [T| was proven in the case when 
the diffusion matrix V"{-) for p,17[) is diagonal. In particular it was shown that 
there is a positive constant C depending only on d, A/A such that 

(1.20) 0< G(z,s;x,t,cb) < -r-, % rr^cxp \x~z\ 

K 1 ~ [A(t - s) + 1]<V 2 [ ^A(t -s) + l 

Hence (|1.19l) . (|1.20l) imply that the uniform intcgrability condition (|1.14[) holds 
for any q with q > 1 + 2/d and the bound on the RHS of (|1.14[) can be taken 
independent of m as m — > 0. Hence if d > 3 the condition (jl . 14[) holds in the limit 
m — > for some q with 1 < q < 2. 

It has been shown by Funaki and Spohn [17] (see also [9]) that if d > 3 then there 
is a unique limit as m — > of the stationary process defined by (|1.15[) , (| 1 . 1 6[) . In 
§3 we shall extend the rate of convergence results in homogenization of (|1.1[) , (|1.2p 
obtained in [11] for the massive field stationary process (|1.15[) , (|1 . 16[) with m > 
to this massless m — > stationary process. In particular we prove the following: 

Theorem 1.1. Let V : R d ->■ R foe a C 2 function such that V"(z), z G R d , is 
a diagonal d x d matrix which satisfies the quadratic form inequality with 
a(-) = V"(-). iei a : R R, d ( d +!)/ 2 & e a (7 1 function on R wii/i wa/ues m the 
space of symmetric d x d matrices which satisfies the quadratic form inequality 
an d has bounded first derivative Da(-) so HDa^Hoo < oo. For d > 3 let 
(f2, J 7 , P) be the probability space of the stationary process </>(•, •) determined by the 
limit as m — > of the stationary process defined by \1.15)) , $1.16]) , and set a(-) 
in ( TO]) to be a(0) = a(^(0, 0)), G fi. Let / : R d -)■ R foe a C°° function of 
compact support, u £ (x,t,Oj) the corresponding solution to U.l\) . M.'0j) with h(x) = 
f(ex), x G Z d , and Uh om (x,t), x G R d ,t > 0, the solution to U.5\) . fl.6\) . Then 
there is a constant a > depending only on d, A/A and a constant C depending 
only on d, A, A, ||-Da(-)||oo, /(•) such that Ji.7| ) holds. 

Remark 1. The exponent a > in \1. ?] ) can foe taken equal to 1 if d > 5 and 
t/ie ratio X/A is sufficiently close to 1. In jllj t/ie matrix V"(-) is not required to 
be diagonal since we use the fact that the Poincare inequality il.lO]) holds for the 
massive field stationary process. In the Gaussian case where V(-) is quadratic Ji.7| ) 
also holds without the restriction that V" {■) be diagonal. This follows from the fact 
that a bound on the Green's function defined by hi. 18^ similar to U.20\) holds in 
this case. Another way of seeing it is to note that the field </>(•, •) is a linear trans- 
lation invariant function of another field a)(-, ■) as in 11.13\) with probability space 
which does satisfy a Poincare inequality. This property of (/)(•,•), being the convo- 
lution of a function with another field whose probability space satisfies a Poincare 
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inequality, does not seem to generalize to the case of uniformly convex V(-) which 
is not quadratic. One reason for this is that the measure for the stationary process 
</>(■, ■) associated with hi. 15}) , hi. lb}) appears to be log concave in </>(•, ■) only in the 
Gaussian case when V(-) is quadratic (see Appendix). In contrast, the invariant 
measure (1.16\) for the stationary process is easily seen to be log concave when V(-) 
is convex. Hence the Brascamp-Lieb inequality [3j implies that a Poincare inequal- 
ity holds for the gradient V </>(■) of the invariant measure field (/)(■) of H.16\) ifV(-) 
is uniformly convex (see |13| ) . 

Parallel to QT| we also establish for the massless field stationary process point- 
wise convergence at large length scales of the averaged Green's function for the 
initial value problem (jl.ip , (|1.2[) to the homogenized Green's function for the initial 
value problem (jl.5|) . (|1.6p . The averaged Green's function G a (x,t), x £ Z d ,t > 0, 
for (HflJ, CL2) is defined by G a (x,t) = ( u(x,t,-) ), where h(-) in (HU) is the 
Kronecker delta function h(x) = if x ^ and h(0) = 1. 

Theorem 1.2. With the same environment as in the statement of Theorem 1.1, 
let G ahoin (x,t), x G R d ,t > 0, be the Green's function for the homogenized prob- 
lem il.5\) . ll.6\) . Then there are constants a, 7 > depending only on d and 
the ratio A/A of the constants A, A of {1-4% and a constant C depending only on 
||-Da(')lloo, A, A, d such that for At > 1, 

(1.21) \G a (x,t)-G^ om (x,t)\ < [At+ ° ](d+a)/2 exp 
(1.22) 

\VG a (x,t)-VG &ho Jx,t)\ < [At + 1] g +1+a)/2 exp 
(1.23) 

\VVG a (x,t)-VVG ahom (x,t)\ < [At+1] ^ +2+a)/2 exp 

The limit as m — > of the invariant measure (|1.16l) is a probability measure on 
gradient fields oj : Z d — >• R d , where formally lo(x) = V<fi(x), x £ Z d . This massless 
field theory measure is crgodic with respect to translation operators [17] for all 
d > 1. In the case d = 1 it has a simple structure since then the variables lo(x), x £ 
Z, arc i.i.d. Note that in the probability space (H,T,P) for the massless field 
theory, the Borel algebra T is generated by the intersection of finite dimensional 
rectangles and the hyperplanes imposing the gradient constraints for w(-). For 
d > 3 the gradient field theory measure induces a measure on fields (j> : Z d —> R 
which is simply the limit of the measures (|1.16j) as m — > 0. For d = 1, 2 the m —> 
limit of the measures (|1.16[) on fields <j> : Z d — > R docs not exist. Naddaf and 
Spencer showed in [31] that the 2-point correlation function for the massless field 
can be represented via the Hclffcr-Sjostrand formula [22] as the expectation value 
of a Green's function for a divergence form PDE with random coefficients. Using 
this fact and techniques of homogenization theory they were able to prove that 
averages of the function x — > ( </>(x)</>(0) ) over large length scales converge to the 
solution of a constant coefficient elliptic PDE 

(1.24) V*a hom V Uhom (a-) = f(x), x e K d . 

Using the techniques of the current paper wc obtain a convergence theorem for the 
correlation function ( <p(x)<p(0) ) which is pointwise in x: 



■ Jl I x 1 



• Ji 1 N 1 
_ 7min | W) __| , 



-7mm < \x\, 
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Theorem 1.3. Let V : R d ->■ R be a C 3 function such that V"(z), z e R d , is 
a diagonal d x d matrix which satisfies the quadratic form inequality {l-4\ l with 
a (') = V"{-) and also ||V'"(-)||oo < °o. Let G ahom (a;), x £ R d , he the Green's 
function for the Naddaf- Spencer PDE \1.2J$ 7 and {■) denote the massless field theory 
measure for the m — > limit of H1.16\) . Then for d > 2 there is a constant a > 
depending only on d and the ratio A/A of the constants A, A of lll.ty) . and a constant 
C depending only on \\V"'(-)\\oo, A, A, d such that for x G Z d — {0}, 

(1-25) | ( <Kx)<j>(0) >-G ahom (x) | < C/\x\ d - 2+a , 

(1.26) | ( V<j>(x)<j>(p) ) - VG ahom (x) | < , 

(1.27) | ( VV*(t>(x)(f>(0) )-VV*G ahom (x) | < C/\x\ d+a . 

Remark 2. Note that both terms ( (j)(x)(j)(Q) ) and G ahoin (x) on the LHS of U.25\) 
are divergent in dimension d = 2. However the difference suitably defined is finite. 
The exponent a > can be taken arbitrarily close to 1 by choosing A/ A sufficiently 
close to 1. 

In the case when V(z) = \z\ 2 /2 + aJ2j=i cosz,,-, z £ R d , the m —> prob- 
ability measure (|1.16j) describes the dual representation of a gas of lattice dipolcs 
with activity a (see [5]). The inequality (|1.25|) was proven by Dimock and Hurd 
(Theorem 2 of [15]) for a continuous space version of the dipole gas under the as- 
sumption that the activity a is sufficiently small. It is not possible to obtain from 
the renormalization group method they use a reasonable estimate on the value of 
a for which their theorem holds. 



2. Variance Estimate on the Solution to a PDE on SI 

We recall some definitions from [TT]. For £ 6 R d and 1 < j < d we define the £ 
derivative of a measurable function ip : SI — > C in the j direction by 9 3 -£, and its 
adjoint by 9*^, where 

(2.1) d 3 ^(u) = e- ie i-«V(Te 3 .,o w)-^(w), 

c^VM = ^■^(r-e,ow)^(^). 

We also define a d dimensional column £ gradient operator 9^ by c\ = (<9i,£, 9^), 
which has adjoint 9| given by the row operator 9| = (9^, ....,9*| A The time de- 
rivative of "0 is defined by 

(2.2) 9V>(w) = lim[V(7D,jw)-V(w)]/<J • 

<5— s-0 

Let %(f2) be the Hilbert space of measurable functions $ : fi — > C d with norm 
ll*llw(n) given by ||\&|[^,™ = ( I^KOli )> where | ■ | 2 is the Euclidean norm on C d . 
Then there is a unique row vector solution $(£, 77, w) = ($i(£, 77, w), 77, w)) 
to the equation 

(2.3) [77 + 9]$(£, 77, u) + 9|a( w )9 c $(£, 77, w) = -0|a(«), 77 > 0, £ € R d , cj e SI, 

such that $(£, 77, -)v £ L 2 (f2) for any v e C d . Furthermore $(£,77, •)« € £ 2 (fi) 
satisfies the inequality 

(2.4) r?||$(£,r?,->||! 2(n) +A]|9^(£,r ?) >|| 2 , w < A>| 2 /A . 
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Letting V denote the projection orthogonal to the constant function, our general- 
ization of Proposition 6.3 of [11] is as follows: 

Proposition 2.1. Suppose a(-) in \2. 3]) is given by &(uj) = a(w(0,0)) where a : 
R n — > Ji d ( d + 1 )/ 2 is a C 1 d x d symmetric matrix valued function satisfying the 
quadratic form inequality {1.4]) and \\Da(-)\\oo < oo. The random field u> : Z^xR^ 
R n is a translation invariant function of a random field uj : Z d x R — > R fc which 
satisfies the uniform integrability condition {1.14]) f or some q with 1 < q < 2. The 
probability space (SI, JF, P) of fields uj : Z d x R — > R fc is assumed to satisfy the 
Poincare inequality {1.10\) . Then there exists qo < 2 depending only on d, A/ A such 
that ifq <q<2 and g G L p (Z d x R, C d ® C d ) with p = 2q/(3q - 2) the inequality 
(2.5) 



xEZ d 



dt g(x, t)d£$(£, ij, t x _ t -)«|| W (n) 



< 



A 



-K u Jg\ 



v G C d , 



ZioZds for a constant C depending only on d, n, k, A/A, qo. 
Proof. From (|1 . 10[) we have that 



/oo 
dt g(x,t)d^, V ,T x ._ t -)v\\ 2 n 
-oo 



xez d 



< 



ds 



z<£Z d ' 



du)(z, 



E 



(n) 



dt g(x,t)d^,r],' 



■)v 



From the chain rule we see that 
(2.7) 

d 

-Q £ *(^7?,T a ,,_f)« = E 



dQ(z, ■ 



dt' 



du)(y,t 



7 rfl c $(^i7,r a , I _f)t 



dtau;(z, -s;y,t',£o) 



Next we do a translation of the functions on the RHS of (|2.7|) . Translation of a 
function G : ft — > C through (x, t) G Z d x R is defined by T x ,tG(ui) — G(t x ^uj), Co G 
SI. For a function G : SI — > C there are two possible notions of translation through 
(x, t), the first being given by t x ^G{uj) = G(r Xt tUi), uj G SI. Since uj is a function 
of a) G ^ we can also define translation through (x, t) by regarding G : SI — > C as 
a function on SI and doing the translation on SI. It follows from the translation 
invariance property (|1.12[) that both of these notions are the same. Now using the 
translation invariance of the probability measure P on S7 we conclude from 

mm that 

dt g(x,t)dt$(£,ri,T x ,-r)v\\^ a) < Kp £ / ds 

-OO . — n fi J -OO 



xez d 



ez d ' 



/oo ,.00 r q 
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We define a function u : Z d x R x O — > C k by 
(2.9) 

/OO 
d*' t'; £, r?, t z ._ s uj)v\ <fou;(0, 0;y + z,t' - s, u) , 

sea- -°° 

where £, 77, w)w : Z d x R -) C™ is the gradient of 77, with respect to 

uj G f2. Observe now from (|1.12p that 



(2.10) d^uj(z, -s;y,t',ti) 







duj(z, — s) 



u(y,t',u)) 







duj(z, —s) 
Hence we have that 
(2.11) 

/OO 
dt' 
-00 



u)(y - z, t' + s, t z - s ui) = dzu(0,0;y - z,t' + s,t z - s Q) . 



y ez d ■ 
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T- Z:S dcjUj(z,-s;y,t' ,Cu) 



J(z-x).£ 



J(z-x).£ 



/OO 
dt' [d^(y - x,t' + t;£,r),T x - z , s - t ui)v] dQ(j(0, 0; y-z, t'+s, ui) = 

IE*" -°° 

/OO 
dt' [d w $(y, t'; £, -q, T x - z , s - t ui)v] dcjUj(0, 0; y + x - z, t! + s - t, ui) 
-00 



u(x — z, t — s, uj) 



It follows from (|2.1ip that (|2.8j) can be rewritten as 

/OO 
dt g(x,t)d^(tV,r x ,-f)v\\ 2 mn) < 
. ., . -OO 

/OO poo 
ds(\ / dtg(x,t)e l{x - z)< Vu(x-z,t~s,-)\ 2 2 ) 

-00 J -co 



z£Z d 



x£Z d 



In [TT| we defined the £ derivative of a measurable function i/j : Z d x R x f! - >C 
in the j direction by -Dj,£, and its adjoint by D* ^, where 

(2.13) Dj t £ip(x,t;oj) = e~ %ei ^^(x - ej,t;T ej uS) - ip(x,t; uj), 

D*^ip(x,uj) = e tej ^ip(x + ej,t; t^ Bj uj) - ip(x,t; uj). 

The corresponding d dimensional column £ gradient operator D$ is then given by 
= (Dx^, Dd,£), and it has adjoint _D| given by the row operator Z?| = 
(£)* ^ , ^)- We also defined the time derivative D of ip : Z d x R x fl — s- C by 

(2.14) D ip(x,t;uj) = \im[ip(x,t - 5;t s uj) - v(j(x,t;uj)}/5 . 

(5— >0 



We see from (|2.1[) . ()2.2j) that these operators satisfy the identities 
(2.15) 



du)(y,t) 
d 



duj(y,t) 



d 6 *P(oj) = D s d^(y,t;uj), y G Z d , t 6 R,w e fl, 
dip((j) = Dodutpiy, t; uj), y € Z d , i e R, uj e O, 
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for differentiable functions ip : fi — > C. A similar relationship holds for the adjoints 
Hence on taking the gradient of equation (|2.3[) with respect to oj(-) we 
conclude from (|2 . 1 5[) that 

(2.16) [ ?? + D ] d u $(y, t'; £, r?, u)v + D|a(w(0, 0))D € d^fo, t'; £, r?, w)« 

= -Dl[6{y,t')Da{u;(0,0)){v + d^^,r],uj)v}} for y G Z d , t' G R, w e O. 

Evidently f|2 . 16[) holds with w G £1 replaced by t z .- s oj for any z e Z d , seR. We 
now multiply (|2.16j) with t z ^ s ui in place of u on the right by e~ lz ^cfc,u;(0, 0; y + 
— s, £>), sum with respect to y G Z d and integrate with respect to t' G R. It 
then follows from (gSJ, (gUp that 

(2.17) 7? tt(z, s, &) - S ' ^ + V*a(u(z, -s))Wu(z, s, u) = -V* f(z,s,u) , 

os 

where the function / : Z d x R x Ct — > C d ® C fe is given by the formula 

(2.18) f(z,s,S>) = Da(u}(z,-s)){v+d^,r],T z ,- s u)v}e- iz <dM^O\z,-s,u}) . 

For any 1 < q < oo we consider the function / as a mapping / : Z d x R — > 
i 2 (Sl, C d <g> C fc ) with norm defined by 



where \\f(y,t',-)\\ 2 is the norm of f(y,f-) G L 2 {S},C d ® C fc ). Now from ([211) it 
follows that d$®(£,rj,-)v € W(Q) and ||<9f$(£, T), -) u ll«(n) < A|v|/A. Hence if the 
inequality (fTTIf holds then the function / is in L q (Z d x R,L 2 (Cl,C d ® C fc )) and 
11/11, < ||-Da(-)||oo(l + A/A)H# w>g . We see from ([2~T7|) that in the case q = 2 then 
V« is also in l?{Z d x R,L 2 (Sl,C d <g) C fe )) and ||Vu|| 2 < H/lh/A. It follows now 
from ([2~T2"]) that (|23|) holds with g = 2 and p = 1. 

To prove the inequality for some p > 1 we use the parabolic version of Meyer's 
theorem [29]. We note that just as the Calderon-Zygmund theorem applies to 
functions with range in a Hilbert space (36], Jone's theorem for parabolic multipliers 
[23] also applies to functions with range in a Hilbert space (see [H]). We conclude 
that there exists qo depending only on d, A/A with 1 < qo < 2 such that if ||/|| q < 00 
for any q satisfying qo < q < 2 then ||Vu|| 5 < 2||/|| 9 /A. Assume now that (|1.14[) 
holds for some q in the interval qo < q < 2. Then by Young's inequality for 
convolutions we see from (|2.12j) that (|2.5j) holds with p = 2q/(3q — 2). □ 



The basic approach of [IT] is to use the fact that the solution to (|1.1|) can be 
expressed by a Fourier inversion formula. For r\ G C denote its real part by K77 £ R 
and its imaginary part by 3?7 G R so that rj — + and similarly denote the 
real and imaginary parts of £ G C d by S£ G R d whence £ = 5R£ + We 
consider solutions to the equation 



[ri + d]&{€,n,w)+'Pd£a{w)dt&{€,Ti,w) = -V%a.(u), 3?77 > 0, £ G R d , u G fi. 



(2.19) 




3. Proof of Theorem 1.1 and Theorem 1.2 



(3.1) 
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As with (|2.3j) there exists a unique solution to (|3.1|) such that $(£,77, ') v £ £ 2 (fi) 
for any v £ C d . Furthermore $(£,77, -)w £ L 2 (f2) satisfies the inequality 

(3.2) 77, >||| a(n) + A||a^(f , r?, -)^II^ C «) ^ ^MV* . 

If £ = the solution $(£,77,0;) to ()2.3f) has zero mean so ( $(0,77, •) ) = 0. Hence 
the solutions to (|2.3[) . (|3.ip coincide if £ = but are in general different. For £ £ R d 
and 77 £ C with 5R?7 > let e(£) 6 C d be the vector e(£) = <%1 and 77) be the 
d x d matrix 

(3-3) g&T?) = (a(-) > + < a(-)0f *(£,»?, •) ) , 

where <&(£, 77,0;) is the solution to (|3.ip . The solution to (11.11) . (]1.2[) is shown in [TU] 

to be given by the formula 

(3.4) 

If the environment (fi,.F, P) is ergodic then the limit lim ?? _ i .o ?(0, 77) = ah om ex- 
ists, and ahom is the diffusion matrix for the homogenized equation (|1.5[) . Let 
^(£,77), £ £ [— 7r,7r] d , 5R?7 > 0, be the Fourier-Laplace transform of the averaged 
Green's function G a (x,t), x eZ d ,t>0, for (fTTjl , (|T2"|) defined by 

OO 

(3.5) G a (e,77) = / dt J2 G a (x,t)eBcp[ia:.f-7«] . 

J ° xez^ 

It follows from Q3.4p that G SL (£,,r]) is given by the formula 

(3.6) G a (£, 77) = l/[77 + e(0*g(f, »?)e(0] for £ G [-tt, 7r] d , ^77 > 0. 

In [TT] it was shown (see especially §3 and Theorem 4.2 of [TT]) that Theorem 1.1 
and Theorem 1.2 are consequences of the following: 

Hypothesis 3.1. For £ £ C d , 77 £ C f/iere exisi positive constants C\ and a < 1 
depending only on d and A/A, suc/i i/ie function <?(£, 77), £ £ R d , 5ft?7 > 0, has an 
analytic continuation to the region 3£| < C\yJWq/k, < 3?77 < A, and 

(3-7) IMS', */)-<?(£, 7,)|| <GA[ \e-t\ a + \(v'-v)/M a/2 

< 5R77 < 5R77' < A, £',£ £ C d with |3£|, |3f£'| < Ciy/ftq/A , 

where C is a constant depending on the environment and the function a(-). 

Here we shall prove that Hypothesis 3.1 holds for the masslcss field theory envi- 
ronment (fi, J 7 , P) of Theorem 1.1. To do this we recall some operators defined in 
[TT] , For any g £ H(£l), let ip(t;,r],uj) be the solution to the equation 

(3.8) j[v + d]il>(t,v,u>) + dZdt1>{t,v,w) = dZg(u>), K77 > 0, £ £ R d , w £ n. 

The operator on is defined by T^ :V g(-) = d^ip(^,r], •). Let G(x,t), x £ 

Z d , /; > 0, be the solution to the initial value problem 

(3.9) 9G(x,t) +v . = ^ ieZ d ii>0i 

at 

G(x,0) = 8(x), x£Z d . 
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It is well known that there exist positive constants C, 7 depending only on d such 
that G satisfies the inequality 



(3.10) G(x,t) + (t + l) 1/2 |VG(a;,t)| + (t + l)\W*G(x,t)\ 



< 



C 



exp 



-7 mm < \x\, 



[t + l] d / 2 

The operator is also given by the formula 



t+1 



for x G Z d , t>0. 



(3.11) T in g(uj) = A e - rit dt ^ {W*G(x,At)}* exp[-ix4] g{r x _ t w) 



< 



It easily follows from p. 81) that T^. l; is a bounded operator on W(£l) with || 2^,r? || -h(o) 
1 provided £ G R d ,5R?/ > 0. Furthermore by Lemma 2.1 of [TT] the function 
(£i ? 7) ~^ horn R d x R to the Banach space of bounded linear operators on 
W(n) has an analytic continuation to a strip < $lr) < A, < C^ftffrj/A where 
C is a constant depending only on d. 

Let b be the d x d matrix valued function b(cj) = Id — a(w)/A, u£fl, whence 
(|1.4[) implies the quadratic form inequality < b(-) < (1 — A/ ' A)Id- It is easy to 
sec from (13.31) that 



(3.12) 



qfor,) = (a(-))-A^(b(.)[PT e >(.)r 



m— 1 



We consider £ G C d ,rj G C with £ having fixed imaginary part, 77 having fixed 
positive real part, and satisfying the conditions of Hypothesis 3.1. For k — 1,2,.., 
we define an operator T^c^s^ from functions g : Z d x R — > C d ® C d to periodic 
functions T k ^ ^ v g : [-7T, Tr] d x R x ->• C rf ® C d by 



(3.13) Tfc,Q«,!R>7 ^ •) 



E 



dt 5 (MK- t Pb(-) [PT^M-)] 



fc-i 



where £ = 77 = afy+^in (|3T3)) . For 1 < p < 00 let LP(Z d xR, C d ®C d ) 

be the Banach space of d x d matrix valued functions g : Z d x R — >• C d ® C d with 
norm ||(7||p defined by 



(3.14) 



\g\\p= sup V / 



dt \g(x,t)v\l 



where \g(x,t)v\2 is the Euclidean norm of the vector g(x,t)v G C d . We similarly 
define the space L°°([— ir, ir] d x R x 51, C d ® C d ) of d x d matrix valued functions 
<y : [— 7r, tt]'' x R x £1 -> C d <g> C d with norm \\g\\oo defined by 



(3.15) 



Iflllc 



sup 

v£C d :\v\. 



sup \\g(C,9,-)v\\ H(n) 

Ce[-7T,7r] d ,6lGR 



Since ||T e < 1 if £ G R d , SR77 > it follows from (I3T31) . (l3~T4l) that if 3£ = 

then T k .^ Mri is a bounded operator from L 1 (Z d x R, C d ®C d ) to L°°([-7r, ir] d x Rx 
f2, C d (8) C d ) with norm HX^Q^sjfrflKoo < (1 — A/A) fc . In the next section we show 
that T k ^,x v is a bounded operator from LP(Z d x R, C d ® C d ) to L°°([-7r, 7r] d x 
R x 0, C d <E) C d ) for some p > 1 in the case of the environment of Theorem 1.1 and 
estimate its norm llT^c^sRjyl^oo- In particular we prove: 
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Lemma 3.1. Let (Q, J-, P) be an environment of massless fields <f> : Z d x R — > R 
with d > 3, and a : R -> R<*(d+i)/2 be as in the statement of Theorem 1.1. Set 
a(0) = a(0(0, 0)), (j) e tt. Then there exists p (A/A) with 1 < p (A/A) < 2 
depending only on d and A/A, and positive constants Ci(A/A), C2(A/A) depending 
only on d and A/A smc/i £/ia£ /or < SJfoy < A, |3f£| < Ci(A/\)yJ$lr)/A, 
(3.16) 

r^iuoc < Ca(A ^*!.y )llo ° ft - a/a)^- 1 )/ 2 [i + c7 2 |^!7(^/a) r 1 

provided 1 < p < po(A/A). 

To complete this section we show how Lemma 3.1 implies that Hypothesis 3.1 
holds. 

Proof of Hypothesis 3.1. We assume that (£,??) and (£',7?') are as m the statement 
of Hypothesis 3.1. Let g : Z d x R ->■ C d ® C d be the function defined by 

(3.17) ff (a:,t) = A{VV*G(a;,Ai)}* e - i ^'-' ? ' t -A{VV*(?(a;,Ai)}* e"^-"* , 

where the Green's function G(; •) is defined by ([321) • It follows from (pTTTj) . (|3T2)) 
and Lemma 2.1 of [TT] that the constant Ci > in (|3.7[) can be chosen depending 
only on d and A/ A so that 
(3.18) 

OO 

\\W,v')-q(Z,ri)}v\\ < CzAj^m^^gm^VrMnn) for |9£|, \$S?\ <C iy /®ii/£, 

k=l 

where C2 is a constant depending only on d, A/ A. We can see from (|3.10[) that there 
is a constant C\ depending only on d such that if < Ci^JWqjA then the 

function g is in L p (Z d x R, C d €5 C d ) for any p > 1. Furthermore if < a < 1 and 
p > (d + 2)/(d + 2 — a) then satisfies the inequality 

(3.19) Hffllp < CpA 1 " 1 ^ |£' - £|» + \(V' V)/M a/2 ] , 

where the constant C p depends only on d, p. The Holder continuity (|3.7[) for suffi- 
ciently small a > follows from (|3.18|) . (|3.19p and Lemma 3.1. □ 



4. Proof of Lemma 3.1 

In |llj we proved that the operator Tj^cjt sjf„ of (|3.13[) is for some p in the range 
1 < P < Po(A/A) a bounded operator from LP(Z d x R, C d ® C d ) to L°°([—k, ir] d x 
R x O, C d ® C d ) if the environment (O, J 7 , P) is the stationary process for the SDE 
(|1.15|) with m > 0. Here we take an alternative approach to proving this result 
which will allow us to study the m — > limit of HT^q^sr^I^oo. We first establish 
an inequality for periodic fields : Q — > R on cubes Q C Z d , and then show that 
we can let Q — > Z d since our estimates are independent of Q. Let L be an even 
integer and Q = Ql denote the lattice points of Z d contained in the cube of length 
L centered at the origin. In the following we identify all points x, y 6 Q with 
x — y = Lek for some k, 1 < k < d. 

As in [TTj the Malliavin calculus (7J [33] is the main tool we use to prove Lemma 
3.1. We assume that V : R d -> R is a C 2 uniformly convex function such that 
a(-) = V"(-) satisfies (|1.4[) and m > 0. Letting B(x,-), x e Q, be independent 
copies of Brownian motion, then the SDE initial value problem 
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(4.1) d<t>(x,t) = — — - ^{V(V^x',t)) + m 2 ^x',t) 2 /2} dt + dB(x,t) 

for x £ Q, t > 0, with <f>(x, 0) = for x £ Q, 

has a unique periodic solution <j)(x,t), x £ Q,t > 0, which is continuous in t > 
with probability 1. Wc denote the function 4> corresponding to a particular realiza- 
tion ui of the white noise process dB(-, ■) as </>(&). Let (f2<g,Mab ^"Q.Mab -fg.Mai) be 
the Malliavin probability space associated with the Brownian motions B(x, •), x £ 
Q. We denote the Malliavin derivative of a function G : Qq Mai ^ C at a point 
S S^Q.Mai by _DMaiG(x, t; w), ieQ,(>0. It is well known (see [7] Theorem 5.4) 
that the Poincare inequality (| 1 . 10[) holds for (fiQ,Mal) -^Q.Mal) -fQ.Mai) with constant 
Kp = 1. Thus we have that 



(4-2) Var[G(-)] < ( \\D Mal G(-;Q)\\ 2 ) n 



Q.Mal 



where || • H2 is the Euclidean norm in L 2 (Q x R + ). 

Let <j> : Z d x R — > R be continuous and consider the terminal value problem for 
the backwards in time parabolic PDE 

(4.3) = ±V*V"(V<f>(y,s))V U (y >S ), s<t,y£Z d , 
u(y,t) = uo(y), y£Z d , 

with solution 

(4.4) u(y,s) = ^2 G(y,s;x,t,4>)uo(x) , s<t,y£Z d . 

It is easy to see that if u (-) £ L 2 (Z d ) then u(-,s) £ L 2 (Z d ) for s < t and 
s ) lli 2 (z d ) fs ll u o( , )llL 2 (z d )- The function G satisfies 

(4.5) y, s ; x ' *> ft = 1 for f G z<i ' E s ; x ' = 1 for x e zd < 

and is non-negative if V"(-) is diagonal. In that case the function x — > G(y, s; x,t,<f>), x £ 
Z d , is the pdf for the position at time t of a continuous time random walk started 
at y at time s. If (/) : Q x R — s- R is periodic we can extend it to a periodic function 
<f> : Z d x R — > R. Let uo : Q — > R be periodic and extend it to a periodic function 
uq : Z d — > R. Then the solution to the periodic terminal value problem (|4.3p is 
given by 

(4.6) u(y,s) = Y G Q {y,s\x,t,4>)u {x) , s<t,y£Q 1 
where Gq is the periodic Green's function 

(4.7) G Q (y,s;x,t) = ^ G (y,s;x + Ln,t,<f>) , s<t,x,y£Z d . 

n£Z d 

It was shown in [TT] that the Malliavin derivative of <f>(x,t,ui), x £ Q,t > 0, is 
given by the formula 

{4.8)D Mai <l>(y,s;x,t,uj) = e- m2 ^/ 2 G Q (y,s;x,t,^(Cu)) for < s < t, 
D Ma ,i<fr(y,s;x,t,uj) = for s > t, y£Q. 
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The solution to (|4.3[) can be written in a perturbation expansion by setting 
V"{z) = A[I d - by(z)], z £ R d , where < by(-) < (1 - X/A)I d in the quadratic 
form sense. Then 

oo 

(4.9) u(y,s) = ^2u n (y,s), yeZ d ,s<t, 

n=0 

where uo(y,s) is the solution to the terminal value problem 

duo^sl = A v , s <t,yeZ d , 
os 2 
u (y,t) = u (y), yeZ d , 

and the u n (y, s), n = 1,2, .., solutions to the terminal value problems 



(4.11) 

du n (y,s) A 
ds 2 



V*Wu n (y,s)- V*b v (V<f>(y,s))Vu n -i(y,s)\ , s<t,yeZ'\ 

u n (y,t) = 0, yeZ d . 



It follows from ([33]) that 

(4.12) u (t/,s) = G{v-z,A(t-s)/2)u (z), y eZ d ,s<t. 

zez d 

Similarly we have that for n > 1, 
(4.13) 

u n (z,r) = V / ds VG(y,As/2)by(V<?i(^+j/,r+s))Vu n _i(2:+j/,r+s) zeZ d ,r<t. 

If we set i*o hi (|4.3p to be given by uo(y) = S(y — x), y 6 Z d , then the perturbation 
expansion (|4.9j) yields a perturbation expansion for the Green's function, 



(4.14) G(y,s;x,t,4>) = ^G n (y,s;x,t,i 



n=0 



where the G n are multilinear in by of degree n. By choosing uo(y) = X)nez d ^(.V ~ 
x—nL), y G Z d , we obtain a similar perturbation expansion for the periodic Green's 
function 

oo 

(4.15) G Q (y,s;x,t,(f)) = ^G ntQ {y,s;x,t,4>) . 

71=0 

Next we consider the inhomogeneous problem 

(4.16) = iv*V-'W(y, S ))V^, S )-/(y, S ), S eR,i,eZ d , 
lim u(y, s) = 0, y £ Z d . 

s— >+oo 

Let / : Z d x R —> C be a continuous function such that dt ||/(-, i)lli 2 (z d ) < 00 • 
From Duhamel's formula we see that the solution to (|4.16[) is given by 

dtG(y,s;x,t,cj ) )f(x,t) , s £ R, y 6 Z d . 
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We can similarly consider the inhomogeneous periodic problem where / :QxR-> 
C is assumed periodic and we extend it to a periodic function / : Z d x R — > C. 
The solution to (|4.16|) is then periodic and is given by the formula 



(4.18) 



dt G Q (y,s;x,t,4>)f(x,t) , s£R,y£Q. 



Lemma 4.1. Assume g : Z d x R — > C is in L 2 (Z d x R) and define for m > the 
function v : Z d x R — > C by 

poo 

(4.19) v(y,s) = / dt e' m2 ^- s y 2 G(y, S ;x,t^)g(x,t) , s£H,y£Z d . 



x<£Z d 



Then v is also in L 2 (Z d x R) and \\v(; QH l^z^ xR) < 2m 2 ||g(-, •)||l, a (Z<»xR)- Cor- 
responding to the perturbation expansion !{4.14\ ) the function v can be written as a 
sum 

(4.20) v = J2 Vn whcrc ll u "(-,-)lli 2 (Z"*xR) < — J (1- VA) n ||ff(-,-)lli»(Z'«xR) ■ 



Proof. It follows from (|4.16|) that v satisfies 



(4.21) 



d m 
ds 



2 1 



v(y,s) 



-V*V"(Wcf>(y,s))Vv(y,s)-g(y,s). 



Multiplying (|4.21[) by v(y, s), summing over y £ Z d and integrating with respect to 

s in the interval —T < s < T we see that 

(4.22) 

,2 



yEZ d 



^ ]T f ds\v(y,s)\ 2 <\Y, {\v(y,T)\ 2 -\v(y,-T)\ 2 }+M £ T ds v(y, s)g{y, s 



y ez" 



yez d ' 



T 
-T 



Now from (|4.19p we have that for s £ R, 



(4.23) \\v(;s)\\ L 2 {zd) 



< 



^ m2t/2 \\g(;s + t)\\ L2{zd) dt 
i 



< 



dt \\g(;S 



\ 2 L 2 (Z d ) 



1/2 



Since g £ L 2 (Z d x R) it follows that the last expression on the RHS of (|4.23j) 
vanishes as s — > oo, whence limT-s-oo ||i>(-, T 1 )||L 2 (z ti ) = 0. We similarly conclude 
that liniT^oo \\v(-, — J 1 )||l 2 (z <1 ) = 0- Letting T — > oo in (|4.22p and using the Schwarz 
inequality we see that if g £ L 2 (Z d x R) then v is also in L 2 (Z d x R) and their 
norms are related as stated. 

To prove (|4.20|) we observe that 



(4.24) 



and that for n > 1 

d m 2 
ths~~ 



(4.25) 



d m 2 



v n (y,s) 



vo(y,s) 



A 



S7*Vv (y,s) -g(y,s), 



A 
2 



V*Vv n (y, s) - V*b v (V0(i/, s))Vv„_i(y, s) 
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Arguing as in the previous paragraph we see from (|4.24|) that 

rn? A 
(4.26) — ||vo(-,-)IL=(z<*xr) + v\\ Vv o(;-)\\l*(z*xr,c«) < ^-all5(-r)lli»cz'«xa) ■ 



Similarly we have from (|4.25[) that for n > 1, 
(4.27) 

771^ A A 

-^-H u "(-,-)lli2 (Z d xR) + -||Vu„(-,-)||i, 2 (zdxR,Crf) ^ ■j( 1 - A / A ) 2 H V ' u r l -l(-, -)Hi2 (Z ^xR.,Crf) ■ 

The inequality in P~2D|) easily follows from (g^SJ), (|4~27j) . □ 

Remark 3. TTie result of Lemma 4-1 holds if Z d is replaced by a periodic cube Q 
with the Green's function G replaced by the periodic Green's function Gq of f^. 7[ ) 
wifft perturbation expansion \4- 15\ ). 

In pj] we considered vector valued functions F(cf>) of fields : Z d x R — > R 
and defined the field derivative of _F at (f> to be the function dF(-, •; 0) with domain 
Z d x R which satisfies 

/oo 
ds dF(y,s;^h(y,s) 

for all continuous functions ft : Z d x R — > R with compact support. Let b : R — )• 
j^d(d+i)/2 k e a (7 1 function taking values in the symmetric d x d matrices such that 
||b||oo + H-Dblloo < oo and define b(-) in (|3.13l) as a function of fields <f) : Z d x R — > R 
by setting b((j>) = b(<£(0, 0)). It follows from (|3~TT|) . ff3TT3|) that for w e C d one has 

?fc,3f ,iRr, 3(5R$, 3?7, </>)u = F(^) where 

(4.29) 

/ <Kj f cxp[-r/(t fc _i-t )-i(.Tfe_i-x )-e] 0(2:0, to) 
i^^ez^- 7 - 00 J 

fc-i 

[] {WG(ij - a?i-i, Afo - t j _ 1 ))}*Pb(0(s o , -<o)) ^b(0(.Tx, -tO) ■ • -Pb^^n, -t fc _i))« . 
i=i 

In (|4.29j) we have extended the domain of the Green's function G(x, t) defined by 
(j3T9]l for i > to t < by setting G(x,t) = 0, t < 0. The operator J - V is 
now any linear operator taking d x d symmetric matrix valued functions b(-) of 
<j> : Z d x R — > R to a constant matrix which has the property that ||(I — 7 ? )b|| 00 < 
su P0:Z d xR^R l|b(0)||txi- We can see from (|3.10[) that there exists G\ > depending 
only on d such that if (£, 77) satisfies the inequality 



(4.30) 0<3fy<A, I Off I < CiV^/A , 
and .g e L 1 ^ x R, C d <g> C d ) then F(0) e C d is bounded by 

(4.31) TOJI^CA^Ubl&p&j] 



-(fc-i). 



where the constant C depends only on d, k. Furthermore F is differentiable in the 
sense of (|4.28|) and the field derivative is given by the formula 
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(4.32) 
dF(y,s\i 



( roc 

A *~M II E J dtj \ exp[-T)(tk-i-to)-i(x k -i-x y£] g(x ,t ) 



fe-i 

3=1 

S(x -y,to + s)Db{(j)(xo, -t )) VbMxu-ti)) ■ ■ ■ Pb(4>(x k - U -t*_i))« + • • • + 



5(x fe _i - tfc_i + s)h((f)(x , -t )) b(0(xi, -*i)) • • ■ Db(0(x fc _i, -t fc _i))v 



for (y,s) £ Z d x R, where S(x,t) = S(x)S(t) is the product of the Kronecker 
and Dirac delta functions. If g is also in L 2 (Z d x R, C d ® C rf ) then dF(-, ■;<£) £ 
L 2 (Z d x R, C d ) and from (|3~T0l) we see that 



(4.33) \\dF(.,.;4>)\\ a < CA^WDh 



ifc-i 



-(fe-i) 



fllla + ®V llfllli | E /" rfi I W*G(x, Ai)| 5 



1/2 



igZ 11 



where C depends only on d,k. 



Lemma 4.2. Let (f2,F, P m ) with m > be the environment of massive fields 
4> : Z d x R -> R denned 6y iTQ3)) . fT76)) and 3 : Z d x R -> C d <g> C d a continuous 
function of compact support. Then there exists C\ > depending only on d such 
that if (£, 77) Zzes m £/ie region Ij4-30{ l the operator of i3.13\) satisfies the inequality 



(4.34) ( g(3t£,$Sri t >| 



< 



E 

yez d 



ds 



E 



dt e- m2{t - s),2 G{y, s; x, t, </>)dF(x, t; 



where G is the Green's function defined by {4-4% an d dF is the field derivative 

EH- 



Proof. It follows from (|4.31[) that the LHS of (|4.34p is the expectation of a bounded 
function. From Lemma 4.1 and (|4.33[) we see that the RHS is also the expectation 
of a bounded function. To prove (|4.34[) we use the Poincare inequality (|4.2|) and 
the formula (|4.8|) for the Malliavin derivative. Thus let Q C Z d be the periodic 
cube with side of length L and <fr(x, t), x £ Q, t > 0, the solution to the initial value 
problem (|4.1[) . For T > we denote by <\>t the periodic field 4>t ■ Q x R — > R 
defined by <pT(x,t) = <f>(x,T + t), x £ Q,t > —T where <fi is the solution to (|4.ip . 
and <pT(x,t) = 0, x £ Q,t < —T. We extend the held <f>T to a periodic field 
(j> T : Z d x R -> R. From flU} we have that for the function F of (j4~2"9]) if y £ Q 
and s > —T then 
(4.35) 

<ft e- m2 (*- s >/ 2 G Q (2/, s; i, i, <fr)dF Q (x, t\ 4>t) , 
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where <IFq is given in terms of (|4.32[) by 

(4.36) dF Q (x,t\4>) = dF ( x + Ln,t;4>) . 

n£Z d 

It is easy to see from ()4.32|) that ||AFg(-, •; 4>)\\oo < C for some constant independent 
of 4> so the RHS of (|435|) is bounded. 

The invariant measure associated with the Markov process defined by the SDE 
(|4.1[) is given by the formula 



(4.37) 



exp 



■£V(V0(x)) 



(xf/2 



J rf(/)(x)/normalization. 



We denote the probability space for the corresponding stationary process 4>(x, t), x € 
Q, t G R by (fig, J-g, Pq.th) and expectation with respect to the measure i"g. m by 
(•)n Qm - Evidently (f2g, .Fg, Pg i?n )is invariant with respect to the translation op- 
erators r x t : fig — > fig, x € Z d ,i 6 R. Our first goal will be to obtain a version 
of the inequality (|4.34[) for the operator 7fc of (|3.13p when the random envi- 

ronment is given by (fig, J-g, Pg. m ). To do this we use the fact (see Appendix A 
for a proof) that for any N > 1, continuous bounded function / : TL N — > C, and 
(xi,ti), (x N ,t N ) e Q x R, 
(4.38) 

lim ( /( T (xi,ti),....,0 r (a;jv,iiv) ) , = ( /( *i)> ••••> ^fciv, ) }~ 



It follows from f|4T38]> that 
(4.39) 



N 



N 



lim / v* J]>b(0 T ( 5 

T ^°° \ 3 = 1 



3' 



-tj)) v 



Z Q,Mal 



v'UvbMxj-tj)) 

3=1 



where I — V on the LHS of (|4.39p denotes expectation with respect to Mal and 
on the RHS expectation with respect to (-)o Q m . We conclude from (|4.39l) . Fubini's 
theorem and the dominated convergence theorem that 
(4.40) 

lim ( |FW>t)| 2 )o Q , MaI = < \m\ 2 )n Q . m = ( g^^M 



Next we see from the Poincarc inequality (|4.2[) and (|4.35p that 
(4.41) 



( I^t)| 2 >, 



< 



E 

y£Z d 



ds 



—T 



E 

x£Z d 



dt e- m2 ^- s ^ 2 G Q (y, s; x, t, cj> T )dF Q (x, t; cfr) 



We assume that L is sufficiently large so that the support of g is contained in Q x R. 
It is easy to see then that dPg(y, s; 4>) is given by the RHS of (|4.32[) with Z d replaced 
by Q and the function G(x, At) replaced by the corresponding periodic Green's 
function on Q. Hence ||<iFg(-, •, <^)||L 2 (QxR,c d ) IS bounded by the periodic version 
of the RHS of (|4.33| . From this we see that limsupg^ z < 1 ||dFg(-, ■,<f>)\\L 2 (QxR,C< 1 ) < 
oo. Observe now that we can argue as in the proof of (|4.40|) to conclude that if 
G n ,Q, n = 0,1,2,.., denote the terms in the perturbation expansion ()4.15|) and 
N > then 
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(4.42) 



/ p oo poo N 

^( E / ds E / *e- m2 ^/ 2 ^G„,Q(y, S ;x,t,0T)^Q^,<;0T 
\ y ez d xeZ d s n=0 





poo 




poo N 




1 ds 


E 


/ dt e -" i2 (*- s )/ 2 G n,Q(v, s; x, t, 4>)dF Q {x, t; cb) 


\ y£Z dJ 


— OO 


xGZ d 


,s n=0 



It follows from (|4.40|) . (|4.42[) and the periodic version of Lemma 4.1 that 



(4.43) ( 








poo 

/ ds 


/•OO 

E / dte- m2 ^/ 2 G Q (y,s;x,t, 0)dF Q (x, t;<f>) 


\ y ez d 


J —oo 


x£Z d Js 



Finally we let Q ->• Z d in (|4.43j) to obtain (|4.34l) . We denote by (-)n m expectation 
with respect to the stationary process defined by (|1.15j) . (| 1 . 16[) . It was proved in 
[17]) (see also [9]) that for any N > 1, continuous bounded function / : R N — > C, 
and (xi,tx), (x N ,t N ) G Q x R, 
(4.44) 

Q 1 ™ <£ ( /( ^l)*!))--)^-^,*^) ) ) QQm = (f{^( x li t l)i-->^( x N,tN))) am ■ 

Hence we can using ()4.44j) argue as with the T — > oo limit to conclude that (|4.34j) 
holds for any m > 0. □ 

For £ G R d and u : Z d — > C we denote by V^u : Z d — > C d the column vector 
V 5 u(z) = [Vi^it(z), Vj, 5 u(z)], z G Z d , where Vj^u(z) = e~ ie '^u(z + ej) - 
u(z), z G Z d , j = 1, .., d. The column operator Vj has adjoint which is a row 
operator. If / G L 2 (Z d x R, C d ) and £ G R d , 7? G C with 9?7i > there is a unique 
solution m(£,77, •) G L 2 (Z d x R) to the PDE 



(4.45) 



n 



d_ 

Or 



u(£, 77, 2, r) + V| V c u(f , 77, 2, r) = V£/(z, r), z G Z d , r G R. 



Furthermore we have that 

K7 

A 



5ft ?7 

(4.46) - r ||u(^T7,-)lli2 (zdxR) + ||V € u(^77,-)||i 2(Z d xRiC d ) < ||/(-) IIl^z-xR.O") 



We also see similarly to p. lip that Vjm(£, 77, •) = T^ :V f(-) where 
(4.47) 

P OO 

T c .„ /(z, r) = A e"" s V {VV*G(j7, As)}* exp[-i V .f] /(z+77, r+s) , z G Z d , r G R. 

From Lemma 2.1 of [TT] we have that u regarded as a function (£, 77) — > L 2 (Z d x R) 
has an analytic continuation to the region (|4.30[) with C± > depending only on d. 
For (£, 77) in this region there is a constant C2 depending only on d such that 
(4.48) 

§ll«(€»»?,-)lli a (z*xR)+llVc«(€,r;,-)lli a (z-xR,c-) < [l+^l 2 / (Sfy/A) ] ||/(0lli a(z -xR,c*) 
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For (j> : Z d x R — >• R a continuous function we extend the corresponding Green's 
function G(y, s; x, t, </>) denned by (|4.4|l for s < t to s > t by setting G(y, s; x, t, 0) = 
when s > t. For £ G R d ,?7 G C,v G C d , with SR77 > let / 2 (£, 77, z, r, 0), z G 
Z rf , r G R with range in C d be defined by 

(4.49) / 3 (f , 77, z, r, 0) = L>b(0(z, -r))« e m2r/2 G(0, 0; z, -r, 0) . 

It is evident from (|4.3p . (|4.4p that / 2 (£, 77, •, •, 4>) is in L 2 (Z d x R), whence we can de- 
fine the function ^7, 2, r, 0) as the solution to (|4.45[) with /(•, •) = /2(C) 77, •> 0)- 
Let b(0) = b(0(O, 0)) and for k = 2,3, .., set 



fc-i 



(4-50) d c F k (^r,,<l>) = [VT^b(-)\ 

where T^ >v is the operator (|3.11|) . We then inductively define functions fk, Uk for 
k = 3, 4.., by the formula 

(4.51) fr({,Ti,z,r,<l>) = Db(^-r))3 4 F fe _i(£,^T*,_ r 0) e mV / 2 G(0, 0; z, -r, 0) 

+ b(0(z,-r)V 4 M fc _i(^,7i,z,r,(/)) , 

where for fc = 3,4,.., the function Uk(^,T),z,r,(j)) is the solution to (|4.45[) with 
/(')') = fk(£,V: ; -A)- Theu fe (£, 77, z,r, (f>) and /*.(£, 77, z, r, 0) for£ G C d are defined 
by analytic continuation from their values when £ G R d . 



Lemma 4.3. Let G be the Green's function defined by l[4-4\ ) an d f or an V k > 2, v G 
C d , iei <£F 6e i/ie function GT3J|). XTien i/iere exists C\ > depending only on d 
such that for (£, 77) m £/ie region i/ie following identity holds: 



/oo 
ds e -" 2 (^)/2 G(z5 _ r; y> S) T _ Z!r(/) ) d F(y, s; r_*, r 0) = 

/OO 
d< «?(x, t) e m2 (*- r )/ 2 G(0, 0; x-z, r-t, ^Db^x-z, r-t))d s F k (£, 77, r x _,, r _ ( 
-OO 

/oo 
g(x, i)b(0(x - z, r - t))V£i£fc(£, 77, x - z, < - r, 0) . 

Proof. We note that the first term on the RHS of (|4.52|) comes from the sum which 
contains S(xo — y,to + s) on the RHS of ( I4.32|) . in which we make the change of 
variables (y, s) O (x, —t). The remaining part of the RHS of (|4.32j) is the same as 

/oo 
dtg(x,t)b(<P(x-z,r-t)) 

/OO 
ds e-™ 2 ( s+r )/ 2 G(z, -r; y, s, T- z . r <f>)dH k (y - x, s + f; T x - x>r -t<f>) , 

where diJfc is the field derivative of the function 
(4.54) H k {4>) = [PT iiV b{-)] k - 1 v . 

Observe now that 

/oo 
ds e-" i2(s +^ 2 G(z, -r; y, s, r^ r <j>)dH k {y -x,s + t; r x _ z ^ t cb) 
-r 
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poo 

= V / ds e-™ 2s / 2 G(0, 0; y, s, 0)dH k (y -( x -z),s+(t- r); r x _,, r _ t , 

so the expression is just a function of (x — z, i — r, 0). We show by induction that 

for k > 2, 

(4.56) 

/"DO 

V e «jb(e,»7,ar,*,0) = Y\ / dse- raa, /2Gi(o,O; V ,a,0)(iff fc (y-a; J a + t;T !e ,_ t 0). 

Hence the identity (|4.52j) follows from (|4.55|) . 

To prove (|4.56[) we first show that it holds for k = 2. To see this we note from 
(|3.11[) that dH2 is given by the formula 

(4.57) dH 2 ( Vl s;4>) = A {VV*G(y, -As)}* cxp[ ?? s - iy4]Db(^y, s))v . 

Now for fc = 2 follows from ([£47)1 . (|4~4T)|) . (|4~5T)) . To do the induction step 

we use the identity 

(4.58) 

dH k+1 (y, s; 4>) = A {VV*G(y, -As)}* expfas - iy.£\Db{<f>(y, s))c\F fc (£, V, r y , s cf>) 
+A / e-' jr dr {VV*G(z,Ar)}*exp[-iz-0)((t)(z,-r))dH k (y-z,s+r-,T z - r (t)). 

From (|4.58l and the induction hypothesis (|4.56j) we conclude that 

POO 

(4.59) / ds e- m2s/2 G(0,0;y, S ,^)dH k+1 (y - x,s + 1;^^) = 

y & Z* J ° 

f(, v [fk+i(t,V,-,-,<t>)](x,t) = V 5 u k+1 (£,ri,x,t,(l>) . 

□ 

In order to prove Lemma 3.1 we shall need to use a parabolic version of Meyer's 
theorem. 

Lemma 4.4. Let H be a Hilbert space and for 1 < p < oo let L p (Z d x R, H) be the 
space of p integrable functions f : Z d x R — 5- T~L where the V norm of f is defined 
by 

/CO 
ds \\f(y,s)\\ p H . 
-oo 

Then there exists a constant C\ > depending only on d such that if (£, rf) lies in the 
region I4-3Q\ ) the operator T^ j?? of is bounded on L p (Z d x R, H) for 1 < p < oo. 

Furthermore there is a constant C*2 > depending only on d such that the norm 
ofT^ tV acting on L p (Z d x R, Ti) satisfies the inequality \\T^ v \\p < [1 + S(p)](l + 
C2|5£;| 2 /[^ 7 7/A]) where the function 8{-) depends only on d and lim p _ > .2 3(p) = 0. 

Proof. This follows from the argument of Lemma 5.2 and Corollary 5.1 of [TT]. □ 

Proof of Lemma 3.1. We assume g : Z d x R — > C d (g> C d is continuous of compact 
support. Then from Lemma 4.2 and Lemma 4.3 we have that 

(4.61) ( \T k ^ v g(M,$Sr,, >| 2 ) < 2 ^ / dr 

zeZd J-oo 
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J2 / dt g{x,t) e^^^GiO^x-z^-t.^Dhi^x-z^r-t^Fki^r,, 



T~x — z,r — t^ 



x£Z d ' 



/oo 
dt g(x, t)h(4>(x — z,r — t))W^Uk(£,, rj,x — z,t — r,(, 
xez," -°° 

Let (f^J 7 , P) be the massless field stationary process corresponding to the limit 
as ni — y of the massive field stationary processes defined by (|1.15[) . (|1.16p . and 
denote expectation with respect to this measure by (-)n- We have from [17] (see also 
[9]) that if d > 3 then for any N > 1, continuous bounded function / : R w — > C, 
and (o:i,ti), {xn^n) G Z d x R, 
(4.62) 

( /( 0(^1, ii), .-,(/>(x N ,t N ) ) ) J2m = ( /( cj)(xi,t{),....,(t>{x N ,t N ) ) ) n . 
Arguing as in Lemma 4.2 we conclude that if d > 3 then 

(4.63) lim ( iTfe,^ .g(^,377,.)«| 2 ) = ( \U M ^ V g(<$%,$r), -)v\ 2 ) . 

Hence to prove Lemma 3.1 it will be sufficient to obtain an upper bound on the 
RHS of (|4.61| which is independent of m as m —> 0. 

For ^eR^eC with Sir) > the functions F fc (£, r?, k = 2,3,...,<f> G O, 
defined by (|4.50[) satisfy the recurrence equations 

(4.64) + 5^3(^,77, ^+^5^2(^)7,0) = AV 8*^(0, 0))v) , 

[ V + d}F k (^ V ^) + Ad^d e F k (^ v ^) = AP5|[b(^0,0))^F fe _i(e,??^)] if * > 2. 
Then as in (|3.2|) we see that ffc(£, ?y, •) £ L 2 (f2) and 

(4.65) ^||F fc «,i?, Oilmen) + IIW(e,'7»-)llw(n) < (1 - A/A) 2 ^ 1 ) | w | 2 . 

The Fk(£,T],(/)) for £ G C d are defined by analytic continuation from the values of 
Fk(£,,r),<fi) when £ £ R d . From Lemma 2.1 of [TT| we see that regarded as a 
function (£, 77) — > L 2 (f2) has an analytic continuation to the region (|4.30p where C\ 
is a constant depending only on d. For (£, 77) in this region there is a constant C2 
depending only on d such that 
(4.66) 

§ll^^r7,-)||i 2(n) +||9 ? F fe fer7,-)||^ ( n) < (1-A/A) 2 ^ 1 ) [l+C 2 |^| 2 / (S,/A) ] 2{k ~ 1] \v\ 2 

We take H to be the Hilbert space U = L 2 (fl,C d ) and for k = 2,3,.., let 
hf. : Z d x R — ► H be the function 

(4.67) fc fc (.z, r, 0) = Db(^(z, -r))9 e F fe (£, r?, r*,_ r 0) e™ 2r / 2 G(0, 0; z, -r, 0) . 

Then the first term on the RHS of (|4.61[) is the square of the norm in L 2 (Z d x R, %) 
of the convolution g * hk-i- It follows from (|1.20[) and (|4.66|) that ifg>l + 2/rf 
then 
(4.68) 

||ftfc|U«(z*xR,w) < C,A- 1 /«|| J Db(.)||oo(l-A/A)( fc - 1 )[l+C 2 |3e| 2 /(5ft f 7/A)] (fc - 1 V'| , 
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where the constant C q depends only on d, q,X/A. Taking q < 2 we have by Young's 
inequality for convolutions that if for p = 2q/(3q — 2) the function g is in L p (Z d x 
R, C d ® C d ) with norm ([3~14"]) then g * h k is in L 2 (Z d x R, H) and 

(4.69) ||<7* hk\\L'{Z^xR,n) ^ C||.9llLP(Z' i xR,C d (8C d ) l|/ l fc||L9(Z' i xR,W) ) 

for a constant C depending only on d. 

To bound the second term on the RHS of ([4.61)1 we show that the function 
(z,r) -> VjUfc(^,i7,2,r,^) is in L q (Z d x R, H). In fact from (fl~20| . (jL49]l we have 
that 

(4.70) ||/2^,»7 > -,-,-)IIl.(z*kr,«) < CgA-^UDbOIUIwl , 

for a constant C g depending only on g, d. Hence Lemma 4.4 implies that 
(4.71) 

HV^^r?, 0IIl«(z*xr,«) < C ? A- 1 /«|| J Db(-)||oc[l+^)] [l+C 2 |^| 2 /(5ft ?? /A) ] | v | . 

From ([4.51)1 we similarly have that for k > 3, 

(4.72) 

ll/k(6 »7, •, •, OIU^z-xr,*) < C,A-V«||Db(.)||oo(l-A/A)( fc - a ) [l+C 2 |^| 2 /(^/A) ] ( *~ 2) |«| 

i "j "j ■)llL9(Z d xR,-H) i 

where C g depends only on q, d. Hence using Lemma 4.4 we have by induction from 
g7Dl) - (jt72l that for k > 2, 

(4.73) ||/fc«,T ?I . 1 . I .)ll£«(Z*xR,W) < 

C q k[l + S(q)] k - 2 A-^WDh^W^l - A/A)( fc - 2 ) [1 + C a |3£| a /(8bj/A) ] (fc_2) |i;| ■ 

The inequality (|3.16|1 follows now from ([4.73)1 and Lemma 4.4 by taking q sufficiently 
close to 2 so that 1 + S(q) < (1 - A/A) -1 / 2 . □ 



5. Proof of Theorem 1.3 

Wc use an identity relating correlation functions for the massive Euclidean field 
theory measure (|1.16|) and expectations of Green's functions for parabolic PDE 
with random coefficients. Let (f2,.F, P m ) be the environment of massive fields 
<f) : Z d x R — > R as in Lemma 4.2 and a(-) in ([1.1)1 be the function a(<fi) = 
V"(V</>(0, 0)), <j> 6 0. Then the two point correlation function for the invariant 
measure ([1.16)1 is related to the expectation of the Green's function for the PDE 
(HU) by 

POO 

(5.1) ( <f>(x)<f>(0) ) = / e-™ 2t G a (x,t) <ft . 

Jo 

The identity ([5.1)1 is implicit in the work of Naddaf and Spencer [31] but was first 
rigorously proven in [18] (see also [9]). 

To prove Theorem 1.3 we observe that the methods used to prove Theorem 1.2 
for diffusion matrices a(-) of the form a.(<f>) — a(0(O,O)) can also be applied for 
diffusion matrices of the form a(</>) = a(V0(O,O)), where a : R d -> R d /( d + 1 )/ 2 [ s 
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C l , satisfies (fl~4|) . and ||£)a(-)||oo < oo. Instead of the bound (fl~20|) we use the 
inequality 

C 

(5.2) \V x G(z,s;x,t,4>)\ < [A(t _ s) + 1](d+W2 ex P 

Evidently (|5.2[) with /3 = is a consequence of (|1.20[) . This is sufficient to prove 
Theorem 1.3 when d > 3. To prove the theorem for d = 2 we need to use the fact 
that /3 can be chosen strictly positive depending only on A/ A. This follows from 
the Harnack inequality P3|. Furthermore f3 can be chosen arbitrarily close to 1 
provided A/ A is sufficiently close to 1. The result follows by letting m — > in (|5.1I) . 

Appendix A. Diffusion Processes with Convex Potential 

Let W : R fe — » R be a C 2 uniformly convex function such that W"(-) satisfies 
the quadratic form inequality XI \ < ]¥"(■) < Alk for some constants A, A > 0. We 
consider the diffusion process : R + — > R fc which is the solution to the SDE initial 
value problem 

(A.l) d<(>(t) = --VW(<f>(t))dt + dB(t), t>0, 0(0) = 0, 

where B(-) is k dimensional Brownian motion. The invariant measure for the SDE 
is given by 

(A. 2) exp[-W(0)] d^/normalization , G R fc . 

We denote the probability space for the stationary process of functions : R — » 
R fe associated with the SDE (|A.1[) and invariant measure (|A.2[) by (fi, J 7 , P), and 
expectation with respect to (SI, J 7 , P) by (-)q. For T > let 0t : [-T, oo) -> R fe 
be defined by 0t(*) = 0(P + 1), where 0(-) is the solution to (|A.1[) . The stationary 
process measure can be obtained by taking the T — > oo limit of 4>t as follows: 

Lemma A.l. Let f : H Nk — > R be a continuous bounded function. Then 

(A.3) lim ( /(0 T (i 1 ),..,0 T (t 7V )) ) = ( /(0(t 1 ),..,0(i JV )) ) n . 

Proof. The diffusion equation corresponding to the SDE (|A.1[) is given by 

(A.4) = -ivW(0)-V«(0,t) + iA«(^,*) , *>0. 

The solution to (|A.4[) with initial data 

(A.5) w(0,O) = Mo (0), 0GR fc , 

can be written in terms of the Green's function G : R fc x R fc x R+ —> R as 

(A.6) u(<M) = / G{4>,<f>' ,t)u {4>') dcj>' , 0eR A \t>O. 

Now it is clear that for a continuous bounded function / : R fc — > R, 

(A.7) (f(Mh))) = [ G(O,0',T + t 1 )/(0') #' . 

Let (-)w denote expectation with respect to the invariant measure (|A.2[) and 
L 2 v ('R k ) the corresponding space of square integrable functions g : R fe — > R with 
respect to {-)w- Letting [•, -}w denote the inner product on L 2 v (H k ) we see from 
(|A.7j) that for any 8 satisfying < S < T + tx, 



\x — z\ 
y/Ht - «5 + 1 
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(A.8) 

</(M*i))> = [h'Mw whcrc/ 1 (0)=G(O,0,J)exp[^(0)] f exp[-W(<j>')] d<f>' 

and/ 2 (<£) = / G{<t>, T + tx- 5) f {(!>') dtf . 

Since H/2II00 < 11/ II 00 it follows that / 2 € L 2 v (R, k ). We can also easily see that 
for 5 > sufficiently small the function /1 is in L 2 ^ r (R k ). Now we use the fact 
that the operator H = — A + '\JW(4>) ■ V is self adjoint non- negative definitive on 
L 2 v (IL k ) and the constant is an eigenfunction of H with eigenvalue 0. From the 
Brascamp-Lieb inequality [3] the operator H acting on the subspace of L 2 v ('R k ) of 
functions orthogonal to the constant is bounded below by A > 0. Furthermore from 
(|A.8|) we have that 

(A.9) ( f(Mh)) > = [fi,e- H(T+t ^/ 2 f] w , 

whence we conclude that 

(A.10) lim ( f(Mti)) ) = ( fM )w{ f(<l>) )w = ( /(</>) )w ■ 

1 —^00 

We have proven (|A.3|) when N = 1. The identity (|A.3[) for N > 1 can be proven 
similarly. Assuming 0<ti <t 2 <" - <£./v, we have that 

(A.ll) (/(<M*i),~,<M*iv)) ) = [fi,e- H ^ T+tl - s ^ 2 g}w where 



g(4>l) = I G((j)i,(f>2,h-tl)---G((f> N - 1 ,(f)N,t N —t N --i)f((t)-L,...,(f>N)d(t>2---d(l)N 

Letting T — > 00 in (jA.llj) we see as before that (|A.3j ) holds. □ 

Next we wish to obtain a representation of the measure for the probability space 
(f2, J 7 , P) for the stationary process associated with the SDE (|A.1|) and invariant 
measure (|A.2[) . First we consider the Gaussian case so there is a symmetric k x k 
matrix and k dimensional vector b with 

(A.12) W(4>) = ^4>*A4> - b*4> , where XI k <A< AI k . 

The SDE (|A.1|) is explicitly solvable when W(-) is given by (|A.12[) with solution 

(A.13) <j>{t) = f e~ A{t - s)/2 [b/2 ds + dB(s)} , t>0. 

Jo 

It is well known that the measure for the stationary process is Gaussian . We can 
use Lemma Al and (|A.13[) to find formulas for the mean and covariance of (/)(■)■ 
Thus we have that 

(A.14) (4>(t)) n = A~\ covnWi),<t>{h)*] = Tfo-h) = A^e^-^/ 2 . 

The Fourier transform of the covariance is therefore given by T{9) = [9 2 + 

A 2 /4] _1 . Hence the Gaussian measure corresponding to the covariance is formally 
given by the expression 



(A. 15) exp 



1 r°° 


d(f>(t) 




dt 



2 1 

+ -\Acj){t) -b\ 2 dt 



J d(j){t) I normalization 
J ten. I 



Evidently the measure (|A.15|) is log concave. 
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We can obtain a representation of the stationary process measure similar to 
()A.15P for general C 2 uniformly convex functions W(-). To see this we write the 
solution of (|A.4[) , (|A.5|) using the Cameron-Martin formula [23] as 
(A.16) 



u(<M) = E 



exp 



1 



VW(B(s)) ■ dS(s) 



\VW(B(s))\ 2 ds }u (B(t)) B(0) = <f> 



2 Jo ° Jo 

where B{) is k dimensional Brownian motion. We rewrite (|A.16j) using Ito's formula 



(A.17) W(B(t))-W(B(0)) 



[ VW{B(s)) ■ dB(s) + \ [ AW(B(s)) ds . 
Jo 2 Jo 



From ()A.16p . (|A.17|I we see that 
(A.18) 



exp 



~AW(B{s)) + ^\\7W{B{s))\ 2 ds \e- w{B(t))/2 uo{B{t)) | B(0) 



The identity (|A.18|) can be alternatively obtained using the Feynman-Kac repre- 
sentation [2H for the solution to the PDE 



(A.19) 

with initial data 
(A.20) 



= VicPM^ + ^Av^t) , t>0, 



w(0,O) = «o(0), 0eR fc . 
Thus we have from the Feynman-Kac formula that 



(A.21) 



v(<t>,t) = E 



exp 



V(B{s)) ds }v (B{t)) B(0)=cf> 



The formula (| A. 18)) follows now from (|A.21j) using the fact that if u(4>, t) is the 
solution to (|A.4jl . (|A.5j) then the function v(cj>,t) = exp[-W(<j))/2]u(<f),t) is the 
solution to (|A.19j) , (|A.20[) with 

(A.22) V{4>) = ^AW{4>) -\\VW{4>)\ 2 , v Q {4>) = e^[-W{4>)/2)u {4>) ■ 

To obtain the representation for the measure of the stationary process we use 
Lemma Al. Thus from (|A.18j) the LHS of (|A.3|) is given by the formula 



(A.23) (f(Mh),..,MtN)) > 

e W(0)/2 E 



exp 



i rT+t N +T' , 1 

--AVF( J B( S )) + i |Vt¥(i?( S ))| 2 ds 



-W(B(T+t N +T'))/2^ B{ r r + ^ gp, + ^ B(T + ^ B ^ = q 



for any T 1 > 0. Recall now that the Brownian motion measure B(s), s > 0, has 
the representation 



(A.24) exp 



d(f>{s) 



ds 



ds 



[ d<j>(s) I normalization with (f>(0) = 
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Hence on letting T, T' — > oo in (|X23|) and using (|A.24p we see that limiting measure 

denned by (|A.3j) with expectation (-)n has the representation 

(A.25) 



cxp 



d<j>{t) 



dt 



lAW(4>(t)) + \\vw(4>(t))f 



dt 



normalization . 



In the Gaussian case (|A.12[) the representation (|A.25|) is equivalent to (|A.15|) since 
AW(-) is constant and is therefore part of the normalization constant. The measure 
(|A.25j) is log concave when W(-) is quadratic, but it is easy to see that even if W(-) 
is a small perturbation of a quadratic the measure is no longer log concave. 

Acknowledgement: The authors would like to thank Tom Hurd and Tom 
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